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Dispersion  Analysis  of  Charleston,  South  Carolina,  Ocean  Dredged  Material  Disposal 
Site(MPDRP-94-l) 


ISSUE:  The  site  designated  for  the  offshore 
disposal  of  material  dredged  from  the  Qiaries- 
ton,  SC,  harbor  is  located  soudieast  of  the  har- 
bOT  entrance.  Historically,  materials  depos¬ 
ited  in  the  Ocean  Dredged  Material  Disposal 
Site  (ODMDS)  have  been  predominately  fine¬ 
grained  sand  with  some  silt  and  clay.  Be¬ 
cause  the  material  from  a  proposed  inner 
harbor-deepening  project  contained  a  higher 
percentage  of  silt  and  clay,  the  Environmental 
Protection  Agency  (EPA)  conducted  a  video 
survey  to  determine  if  the  ODMDS  contained 
ecologically  sensitive  areas  that  might  be  af¬ 
fected  by  the  disposal  of  the  finer-grained  tiu- 
terial.  The  EPA  identified  areas  of  live  bot¬ 
tom  habitat  within  the  Charleston  ODMDS. 

RESEARCH:  The  dispersive  characteris¬ 
tics  of  the  disposal  site  were  investigated  as  a 
fiinctimi  of  die  local  wave  and  current  environ¬ 
ment  Site-specific  wave  and  current  informa¬ 
tion  was  used  for  the  boundary-condition  de¬ 
velopment  for  numerical  modeling  of  the 
short-term  effects  of  the  disposal  operation  as 
well  as  the  long-term  bdiavior  of  the  disposal 
mound.  Results  of  this  study  contribute  to  the 
body  of  knowledge  from  siinilar  site- 


characterization  studies  and  to  the  credibility 
oof  numerical  modeling  efforts. 

SUMMARY:  Short-term  numerical  simula¬ 
tions  indicated  that  even  when  worst-case  cur¬ 
rents  are  reef  directed,  the  finer-grained 
dredged  material  does  not  appear  to  pose  a 
threat  if  the  disposal  location  is  situated  at 
least  1.5  miles  from  die  bottom  habitat.  Long¬ 
term  dredged  material  ntound  movement  is 
relatively  small  under  normal  conditions. 
However,  storm-induced  currents  may  cause 
erosion  and  transport  of  material  that  could 
pose  a  threat  to  tire  coral  reefs  if  conditions 
are  severe  and  storm  durations  are  significant 
unless  disposal  mounds  are  located  a  reason¬ 
able  distance  to  the  northeast  of  the  area  of 
bottom  habitat 

AVAILABILITY  OF  REPORT:  The  report 
is  available  through  the  Interlibrary  Loan  Ser¬ 
vice  from  the  U.S.  Army  &igineer  Waterways 
Experiment  Station  (W^)  Library,  telephone 
number  (601)  634-2355.  National  Technical 
Information  Service  (NTIS)  report  numbers 
may  be  requested  from  WES  Librarians. 

To  purchase  a  copy  of  the  report,  call  NTIS  at 
(703)487-4780. 
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PREFACE 


The  site  designation  methodology  used  in  this  report  was  developed  at 
the  Coastal  Engineering  Research  Center  (CERC)  and  the  Hydraulics  Laboratory, 
US  An^  Engineer  Vaterways  Experiment  Station  (WES)  by  work  units  of  Technical 
Area  1  (TAl) .  Analysis  of  Dredged  Materials  Disposal  in  Open  Water,  of  the 
Dredging  Research  Program  (DRP)  of  Headquarters,  US  Army  Corps  of  Engineers 
(HQUSACE).  Messrs.  Robert  Casqtbell  and  Glenn  R.  Drummond  were  DRP  Chief  and 
TAl  Technical  Monitor  from  HQUSACE.  respectively.  Mr.  E.  Clark  McNair,  Jr., 
CERC,  was  DRP  Program  Manager  and  Dr.  Lyndell  Z.  Hales,  CERC,  was  Assistant 
Program  Manager.  Dr.  Nicholas  C.  Kraus,  Senior  Scientist,  Research  Division 
(RD) ,  CERC  was  Technical  Manager  for  DRP  TAl . 

This  report  describes  an  investigation  of  the  dispersion  characteristics 
of  the  existing  Charleston  Ocean  Dredged  Material  Disposal  Site,  located 
southeast  of  the  entrance  to  Charleston  Harbor.  The  study  was  conducted  by 
CERC  at  the  request  df  the  US  Army  Engineer  District  (USAED) ,  Charleston.  Two 
sources  of  current  measurements  were  used  in  the  present  analysis.  The  first 
source  was  data  collected  by  the  National  Oceanic  and  Atmospheric  Administra¬ 
tion  (NOAA)  as  a  component  of  the  Charleston  Harbor  Oceanography  Project. 
Appreciation  is  extended  to  Drs.  Henry  Frey  and  Wayne  L.  Wilmot  and  Mr.  Cary 
R.  Wong,  NOAA,  for  providing  this  inforoiation.  The  second  source  of  data  was 
provided  by  the  US  Environmental  Protection  Agency.  Appreciation  is  extended 
to  Mr.  Philip  J.  Murphy,  EPA,  for  making  these  data  available.  Appreciation 
is  also  extended  to  Mr.  Steve  J.  Morrison,  USAED,  Charleston,  for  providing 
Information  and  assistance  crucial  to  the  timely  completion  of  this  project. 

Both  phases  of  the  numerical  investigation  and  the  final  report  were 
prepared  by  Drs.  Norman  W.  Scheffner  and  James  R.  Tallent,  Coastal  Processes 
Branch  (CPB) ,  RD,  CERC.  WES. 

Providing  general  supervision  were  Dr.  James  R.  Houston  and  Mr.  Charles 
C.  Calhoun,  Jr.,  Director  and  Assistant  Director,  respectively,  CERC;  direct 
supervision  of  the  project  was  provided  by  Mr.  H.  Lee  Butler,  Chief,  RD,  and 
Mr.  Bruce  A.  Ebersole,  Chief,  CPB. 


At  th«  tlM  of  publication  of  this  import,  Director  of  WES  was 
Dr.  Robert  V.  Vhalln.  Consiander  was  COL  Bruce  K.  Howard,  EN. 


For  additional  inforaation  on  this  work,  contact  Dr.  Norman 
Scheffner,  (601)  634-3220,  or  Mr.  E.  Clark  McNair,  Jr.,  at 
(601)  634-2070. 
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CONVERSION  FACTORS,  NON-SI  TO  SI  (METRIC) 
UNITS  OF  MEASUREMENT 


Mon- SI  units  of  neasurenent  used  In  this  report  can  be  converted  to  SI 
(oetric)  units  as  follows: 


- WvUiply -  Bv 

cubic  yards  0.7645549 

degrees  (angle)  0.01745329 

feet  0.3048 

niles  (US  nautical)  1.852 

nlles  (US  statute)  1.609347 


— To  Obtain 

cubic  metres 

radians 

metres 

kilometres 

kilometres 
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SUMMARY 


In  this  report  the  dispersion  characteristics  of  dredged  material  place¬ 
ment  operations  at  the  Charleston,  SC,  Ocean  Dredged  Material  Disposal  Site 
(ODMDS)  are  investigated.  The  primary  focus  of  the  study  is  to  determine  if 
material  deposited  at  the  designated  disposal  site  will  migrate  to  the  live 
coral  reefs  that  were  recently  discovered  within  the  boundaries  of  ODMDS. 

This  study  was  conducted  at  the  request  of  the  U.S.  Army  Engineer  District, 
Charleston. 

A  disposal  site  can  be  classified  as  dispersive  or  nondispersive  depend¬ 
ing  on  whether  sediment  is  transported  out  of  or  remains  within  the  designated 
limits  of  the  site.  The  dredged  material  dispersion  characteristics  of  the 
Charleston  ODMDS  were  investigated  in  two  phases,  a  short-term  and  long- tens 
phase.  In  the  short-term  phase  the  potential  effect  of  the  actual  barge  dis¬ 
posal  operation  on  the  local  environment  was  investigated.  This  phase  of  the 
study  represents  the  initial  minutes  to  hours  following  the  disposal  operation 
during  which  time  the  material  is  entrained  and  dispersed  as  it  descends 
through  the  water  column  to  be  deposited  on  the  ocean  floor.  Efforts  were 
focused  on  modeling  the  time  rate  of  change  of  suspended  sediment  concentra¬ 
tion  and  the  total  sediment  deposition  pattern  on  the  ocean  bottom.  In  the 
second  phase,  the  long-term  analysis  focuses  on  the  extent  and  probable  direc¬ 
tion  in  which  local  waves  and  currents  erode  and  transport  the  dredged 
material  mound.  The  methodologies  used  to  accomplish  these  goals  are  thor- 
ou^ly  discussed  in  the  report. 

Short-term  numerical  simulations  were  perfoxnned  for  worst-case  wave  and 
current  conditions.  Results  include  the  water  column  spatial  distribution  of 
the  sand  and  silt -clay  components  of  the  sediment  load  in  the  form  of  sediment 
concentration  (parts  per  million)  above  the  background  level.  Computational 
results  indicate  that  a  significant  fraction  of  the  sand  and  silt/clay 
materials  falls  rapidly  to  the  ocean  floor  and  does  not  impact  regions  beyond 
one -quarter  mile  from  the  point  of  disposal.  However,  a  small  amount  of 
silt/clay  material  remains  in  the  water  coliimn  for  several  hours  after  the 
disposal  operation.  This  cloud  of  suspended  material  is  transported  about  1 
mile  from  the  disposal  point.  The  maximum  thickness  of  the  final  deposition 
was  approximately  0.5  ft,  covering  an  approximate  700-ft-diameter  area.  Depo¬ 
sition  is  confined  to  this  immediate  area. 
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Results  of  the  long-term  simulations  Indicate  that  the  mound  Is  disper¬ 
sive  with  respect  to  normal  wave  and  tldal/clrculatlon  currents,  with 
migration  rates  as  large  as  60  ft/month.  Storm  events  significantly  Increased 
mound  movement.  The  simulation  of  a  moderate- Intensity  event  with  a  24-hr 
duration  showed  the  movind  migration  to  be  approximately  155  ft. 

Based  on  the  findings  of  this  report  It  Is  concluded  that: 

1)  when  the  worst-case  currents  are  reef  directed,  neither  the  sand  nor  the 
sllt/clay  isaterlals  appear  to  pose  a  threat  If  the  disposal  site  Is  situated 
at  least  1.5  miles  from  the  reef,  and  2)  long-term  dredged  material  mound 
movement  Is  relatively  small  under  normal  conditions;  however,  storm- Induced 
movement,  which  can  be  directed  In  virtually  any  direction,  may  pose  a  threat 
to  the  coral  reefs  If  conditions  are  severe  and  storm  durations  are 
significant. 
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PART  I :  INTRODUCTION 


BftgKgrflvmd 

1.  The  designated  site  for  the  offshore  disposal  of  dredged  material 
for  Charleston,  SC,  Is  located  southeast  of  the  jettled  entrance  to  Charleston 
Harbor,  Figure  1.  The  approximately  5*  by  9 -km  designated  site  has  been  In 
use  for  the  disposal  of  dredged  material  since  the  mld-1970's.  Historically, 
materials  deposited  at  the  site  were  predominately  fine-grained  sands  with 
some  silts  and  clays.  The  proposed  Inner  harbor  deepening  project  will 
require  the  disposal  of  approximately  3.0  million  cubic  yards*  of  material  In 
the  existing  Ocean  Dredged  Material  Disposal  Site  (ODNDS) .  Since  the  material 
from  this  project  contains  a  higher  percentage  of  fine-grained  silts  and  clays 
than  the  original  Charleston  Harbor  deepening  project,  an  Investigatory  video 
survey  was  undertaken  by  the  Environmental  Protection  Agency  (EPA)  to  deter¬ 
mine  whether  the  site  contained  ecologically  sensitive  areas  that  may  be 
affected  by  the  disposal  of  fine-grained  material.  Results  of  that  survey 
Identified  "extensive  live  bottom  habitat  areas  within  the  Charleston 
OmDS.”**  The  purpose  of  this  study  Is  to  use  recently  obtained  prototype 
current  data  to  quantify  the  local  current  patterns  and  magnitudes  and  to  use 
these  data  to  Investigate  the  potential  effect  of  the  new  disposal  operation 
on  the  adjacent  environmentally  sensitive  areas. 

abl9£.t.LYe 

2.  The  objective  of  this  study  Is  to  Investigate  the  dispersive  charac¬ 
teristics  of  the  proposed  site  as  a  fvinctlon  of  the  local  wave  and  current 
envlroniBent .  In  this  manner,  the  potential  effect  of  the  disposal  site  on  the 


*  A  table  for  converting  non-SI  units  of  measurement  to  SI  (metric)  units 
Is  presented  on  page  5. 

**  Environmental  Protection  Agency,  14  February  1990,  Memorandum  Concerning 
the  Location  and  Extent  of  Live  Bottom  Areas  In  the  Charleston  Ocean 
Dredged  Material  Disposal  Site. 


Figure  1.  Charleston  ODMDS  location 


underwater  comnunltles  can  be  systematically  Investigated  In  a  two>phase 
numerical -modellng>based  approach.  First,  the  short-term  effects  of  the 
dredging  operation  are  Investigated  to  determine  the  time  rate  of  change  of 
suspended  sediment  concentration  In  the  water  column  as  the  descending  sedi¬ 
ment  plume  disperses  and  Is  transported  from  the  site  by  ambient  currents. 
Additionally,  the  total  sediment  deposition  pattern  Is  computed  to  Indicate 
the  maximum  spatial  extent  of  deposition  on  the  ocean  bottom.  The  modeling  of 
this  short-term  phase  of  the  operation  Is  performed  by  the  Disposal  from  an 
Instantaneous  Dump  (DIFID)  numerical  model  (Johnson  1990) .  This  model  com¬ 
putes  the  convective  descent  and  dynamic  collapse  of  the  sediment  following 
Its  release  from  the  barge. 

3.  The  second  phase  of  the  Investigation  examines  the  behavior  of  the 
sediment  mound  over  long  periods  of  time.  Althou^  It  Is  recognized  that  the 
site  Is  a  dispersive  one,  this  long-term  analysis  focuses  on  the  extent  and 
probable  direction  In  which  local  waves  and  currents  erode  and  transport 
deposited  material  from  the  sediment  mound.  These  simulations  are  performed 
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with  m  eeuplad  hydrodynawlc ,  ••diwant  transport,  and  bathyaatry  changa  aodal 
(Schaffnar  1992  and  in  preparation)  that  conputes  nound  stability  as  a  func¬ 
tion  of  Bound  coapositlon  and  environnental  forcings.  Both  Bodeling  efforts 
require  site-specific  infoniatlon,  including  waves,  currents,  bathyaetry, 
sedlaent  coaposition,  and  disposal  aethods. 

4.  A  realistic  analysis  of  the  dispersion  characteristics  of  the 
Charleston  OMDS  can  be  Bade  only  if  the  coaputations  are  based  on  site- 
specific  wave  and  current  inforaation.  This  investigation  is  fortunate  in 
that  recently  obtained  current  data  for  several  locations  near  and  within  the 
disposal  site  are  available.  Current  aeasureaents ,  collected  by  the  National 
Oceanic  and  Ataospheric  Adalnistration  (NOAA)  as  a  coaponent  of  the  Charleston 
Harbor  Oceanography  Project,  were  provided  to  the  Coastal  Engineering  Research 
Center  (CERC) ,  US  Arnqr  Engineer  Waterways  Experiment  Station  (WES) ,  for  analy¬ 
sis  and  use  in  this  study.  One  current  meter  was  placed  within  the  limits  of 
the  ODNDS  by  the  EPA  to  provide  information  on  local  ctirrents.  These  data 
were  also  supplied  to  CERC  for  use  in  this  study. 

5.  This  report  concentrates  on  the  three  primary  components  of  the 
study:  botindary  condition  development,  short-term  modeling,  and  long-term 
modeling.  The  most  important  component  of  the  three  is  the  development  of 
realistic  boundary  conditions  at  the  site.  The  accuracy  and  credibility  of 
the  numerical  modeling  effort  are  dependent  on  the  realistic  appr*^*- ■  uttion  of 
waves  and  currents  at  the  disposal  site.  The  importance  of  this  aspect  of  the 
study  has  been  stressed  in  similar  site  characterization  studies  and  will  be 
the  subject  of  Part  II  of  this  report. 
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PART  II;  WAVE  AND  CURRENT  BOUNDARY  CONDITIONS 


6.  The  short-term  modeling  phase  of  this  analysis  requires  the  specifi¬ 
cation  of  a  local  velocity  field  representative  of  the  conditions  found  at  the 
disposal  site.  Because  the  DIFID  model  applies  only  to  the  time  ImiBedlately 
following  disposal,  a  single-value  velocity  Is  specified  to  represent  a 
"worst-case  dispersion  scenario"  in  «dilch  the  currents  are  directed  from  the 
disposal  site  to  the  live  bottom  habitat  regions.  In  this  manner,  a  single - 
value,  depth- averaged  velocity  can  be  used  to  quantify  the  potential  effect  of 
the  disposal  operation  on  the  reefs. 

7.  The  long-term  modeling  phase  requires  a  more  nearly  complete  and 
comprehensive  description  of  local  waves  and  currents  because  the  modeling 
approach  Investigates  the  behavior  of  a  disposal  mound  over  long  periods  of 
time  (on  the  order  of  months) .  As  such,  a  realistic  representation  of  the 
wave  and  current  time  series  Is  required  for  the  site;  otherwise  realistic 
predictions  of  mound  stability  cannot  be  made.  The  remaining  paragraphs  of 

II  will  concentrate  on  defining  the  wave  and  current  boundary  conditions 
required  as  Input  to  both  the  short-  and  long-term  models. 

HflYfi.. Height.  Period,  and  Direction  Time  Series 

8.  The  long-term  transport  model  computes  sediment  transport  as  a  fxinc- 
tion  of  time  series  of  both  waves  and  currents.  The  wave  time  series 
coaq>onent  Is  specified  as  a  statistical  simulation  of  the  20-year  hlndcast 
data  base  of  the  Wave  Information  Study  (WIS)  (Jensen  1983),  Phase  III,  Sta¬ 
tion  117  "sea"  conditions.  The  location  of  Station  117  Is  shown  In  Figure  2. 
The  statistical  approach  to  defining  time  series  of  wave  height,  period,  and 
direction  for  a  specific  WIS  station  Is  reported  in  detail  by  Borgman  and 
Scheffner  (1991) .  The  approach  allows  the  user  to  simulate  wave  sequences 
that  preserve  the  statistical  qualities  of  the  entire  20-year  data  base. 
Including  seasonality  and  wave  sequencing.  The  statistically  based  time 
series  provides  a  site  specific  wave  climate  that  is  Ideal  for  the  long-term 
simulation. 

9.  A  1-year  time  series  of  waves  was  generated  as  input  for  the  long¬ 
term  model.  Plots  of  the  simulated  sequence  of  wave  hel^t,  period,  and 
direction  are  shown  in  Figure  3.  To  demonstrate  the  visual  similarity  between 
the  simulated  wave  field  and  actual  hlndcast  data.  Figure  4  represents  a 
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Figure  2.  WIS  hlndcast  data  station  locations 


l*year  tiae  series  of  WIS  data  for  the  year  1956.  Both  Figures  3  and  4  begin 
on  1  Janxiary  and  extend  through  31  Deceaber.  The  siailarity  in  patterns  of 
increased  winter  activity  with  a  decrease  in  intensity  during  the  summer 
months  can  be  seen  in  both  plots.  A  quantitative  comparison  of  the  data  can 
be  seen  in  the  percent  probability  histogram  plots  shown  in  Figure  5,  in  which 
the  probability  statistics  of  the  simulated  wave  height,  period,  and  direction 
are  overlaid  with  those  of  the  1956  VIS  data.  A  compilation  of  computed  maxi¬ 
mum,  minimum,  average,  and  standard  deviation  values  for  both  the  simulated 
data  and  the  hlndcast  data  for  the  years  of  1956,  1960,  1964,  and  1968  are 
shown  in  Table  1.  The  similarity  in  values  demonstrates  that  the  simulated 
time  series  are  statistically  similar  to  the  hlndcast  data  base.  Also 
included  in  Table  1  is  the  direction  from  which  the  greatest  ntmiber  of  waves 
originated.  The  reported  angles  are  with  respect  to  a  shore -parallel  baseline 
so  that  90  deg  indicates  a  wave  field  directed  onshore.  Computations  indicate 
that  the  ujority  of  waves  are  from  the  northeast  and  southwest.  This  obser¬ 
vation  is  consistent  with  the  bimodal  directional  distribution  of  the  VIS  data 
shown  in  the  directional  histogram  of  Figure  5. 
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Tabl*  1 

(kMDTlaon  of  W«v  Statlatica 


Parameter 

Simulated 

1956 

1960 

1964 

1968 

Maximum  wave  height,  m 

3.7 

3.3 

3.3 

3.3 

2.9 

Minimum  wave  height,  m 

0.0 

0.0 

0.0 

0.0 

0.0 

Average  wave  height,  m 

0.5 

0.5 

0.6 

0.6 

0.3 

Standard  deviation,  m 

0.4 

0.4 

0.5 

0.6 

0.4 

Maximum  wave  period,  sec 

10.5 

9.0 

9.0 

9.0 

9.0 

Mlnlmm  wave  period,  sec 

0.0 

0.0 

0.0 

0.0 

0.0 

Average  wave  period,  sec 

3.9 

3.8 

4.4 

4.4 

2.9 

Standard  deviation,  sec 

1.5 

1.8 

1.7 

1.7 

1.9 

Most  common  direction,  deg 

10-20 

20-30 

160-170 

160-170 

0-10 

10.  Station  117  rapresents  a  Phase  III  VIS  hindcast  station,  and  as 
such,  the  hindcast  Is  developed  for  10  m  of  water.  The  following  relation¬ 
ships  were  \ised  to  transform  the  wave  height  from  10  m  to  deep  water  and  then 
to  shoal  the  wave  from  deep  water  back  to  the  disposal  site.  This  transfor¬ 
mation  maintains  the  proper  wave  height  relationship  for  the  wave  shoaling  in 
response  to  the  presence  of  the  disposal  site  (Ebersole,  Clalone,  and  Prater 
1986): 

H  -  Ho  k.  (1) 

where  the  shoaling  coefficient  k.  Is  defined  as: 


The  parameters  k  and  h  represent  the  wave  number  and  local  depth, 
respectively. 

Depth-Averased  Current  Time  Series 

11.  The  analysis  of  the  velocity  patterns  In  the  vicinity  of  the  ODMDS 
begins  with  an  analysis  of  the  multiple -depth  NOAA  data,  since  that  data 
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provides  Inforastlon  on  current  verietlons  with  respect  to  depth.  Results 
will  then  be  supplenented  end  enhanced  by  the  single 'gage  EPA  data  located 
6  ft  off  the  bottom  within  the  boundaries  of  the  disposal  site. 

12.  The  current  information  was  recorded  with  bottom- moored,  self- 
contained,  acoustic  doppler  profilers  manufactured  by  RDI  Instruments  Inc. 

Data  were  provided  for  the  four  locations  shown  in  Figure  6.  Summary  statis¬ 
tics  of  each  instr\uient  location  are  given  in  Table  2. 

13.  Since  depth- averaged  currents  are  used  for  the  present  analysis,  a 
depth -averaged  value  at  each  time-step  was  computed  by  averaging  the  values 
for  the  depth  through  the  water  column  at  each  10 -sec  time -step  of  the  time 
series  shown  in  Table  2.  A  six-constituent  harmonic  analysis  was  performed  on 
the  depth-averaged  data  to  provide  an  indication  of  tidal  influence  in  the 
tidal  signal.  Results  of  the  analysis  are  summarized  in  Table  3.  The 
reported  reduction  in  variance  (RV)  values  indicates  the  degree  of  tidal 


Figure  6.  Station  locations  for  current  meter  data 
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Table  2 

Sunnarv  Statistics  of  NOAA  Current  Data 


Depth 

Mean  Velocity 
cm/sec 

Standard  Deviation 
_ cm/sec _ 

m 

u 

V 

V 

V 

Statign C-2? 

Location  (lat/long) :  79.7923/32.7150 

Samoline  oeriod: 

17  Mav  1988  - 

21  June  1988 

0.9 

18.7 

-6.1 

30.3 

22.1 

1.9 

17.0 

-6.0 

29.1 

21.6 

3.0 

13.7 

-4.1 

26.8 

19.8 

4.1 

10.7 

-3.0 

24.8 

18.7 

5.2 

7.9 

-2.3 

22.8 

18.0 

6.3 

5.0 

-1.8 

20.9 

17.4 

7.4 

2.2 

-1.5 

19.3 

16.7 

8.4 

-0.4 

-1.2 

18.0 

15.8 

Station  C-24 

Location  aat/long> :  79.7260/32.6830 

Samoline  oeriod: 

17  Mav  1988  - 

21  June  1988 

0.3 

7.9 

2.4 

23. 

18.2 

1.3 

6.0 

2.6 

18. 

12.8 

2.3 

6.9 

2.3 

19, 

12.8 

3.3 

6.1 

2.3 

18. 

12.7 

4.3 

4.9 

2.3 

17. 

12.3 

5.3 

3.5 

2.5 

16. 

12.0 

6.3 

2.2 

2.8 

14. 

11.5 

7.3 

1.1 

3.0 

13. 

11.5 

8.3 

0.3 

3.0 

12. 

10.2 

9.3 

0.8 

3.0 

12. 

9.9 

Station  C-33 

Location  (lat/lone):  79.8100/32.7243 

Samoling  Period: 

21  June  1989 

-  18  Seo  1989 

0.0 

23.0 

-5.9 

50. 

37.1 

1.1 

17.5 

-7.0 

48. 

29.1 

2.2 

15.3 

10.9 

49. 

29.1 

3.3 

11.6 

11.4 

48. 

28.0 

4.4 

7.5 

-9.8 

46. 

26.5 

5.5 

3.5 

-9.5 

44. 

24.7 

6.6 

0.1 

-9.5 

41. 

22.8 

7.7 

-2.1 

-9.9 

38. 

20.6 

(Continued) 

Note:  Station  locations  indicated  by  latitude  and  longitude. 
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Table  2  (Concluded) 


Mean  Velocity  Standard  Deviation 


Ma 

26.17 

20.27 

16.31 

7.24 

57.82 

29.59 

18.54 

9.41 

Si 

3.73 

3.15 

1.77 

2.36 

7.73 

4.11 

2.37 

1.41 

N2 

5.09 

4.51 

1.01 

1.72 

9.97 

3.52 

3.22 

1.76 

Ki 

3.70 

0.42 

0.88 

2.15 

5.95 

3.34 

1.88 

0.24 

0i 

5.11 

2.47 

2.05 

1.31 

4.06 

2.00 

1.34 

0.50 

Pi 

8.72 

4.29 

3.10 

1.92 

2.18 

0.85 

1.40 

0.10 

RV 

0.72 

0.73 

0.5 

0.25 

0.9 

0.71 

0.64 

0.35 

Ave. 

9.34 

-3.24 

3.89 

2.61 

9.54 

9.12 

4.30 

2.33 

*  Const.  -  constituent 
ave.  -  average 
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Influance,  a  valua  of  1.0  indicating  that  the  selected  six  constituents  coa- 
pletely  describe  the  time  series.  Conversely,  a  value  of  0.0  indicates  that 
the  tine  series  is  uncorrelated  with  tides.  As  indicated,  the  signals  are 
tide  dominated  and  are  primarily  semidiurnal  in  period. 

14.  A  visual  indication  of  the  nontidal  component  of  the  signal  can  be 
demonstrated  through  the  application  of  a  low-pass  nonrecursive  digital  filter 
that  separates  the  high-frequency  (tidal)  component  of  the  series  from  the 
low-frequency  (storm  or  large-scale  circulation)  component.  Figures  7  through 
10  present  these  data  in  the  form  of  raw  time  series  and  the  U  (east)  and  V 
(north)  component  with  the  superimposed  low-frequency  component.  A  current 
direction  histogram  is  also  shown. 

15.  Since  the  directional  histogram  provides  no  information  on  magni¬ 
tudes,  the  direction  vector  of  each  U,V  pair  was  computed  to  indicate  the 
coupling  of  magnitude  and  directionality.  These  relationships  are  shown  in 
the  form  of  continuous  vector  plots  of  each  time  series  (at  a  60-min  interval) 
and  are  Illustrated  in  Figure  11. 

16.  An  examination  of  the  data  from  Stations  C-23  and  (especially)  C-34 
shows  that  these  current  station  locations  are  dominated  by  the  effects  of  the 
jetty  entrance  and  therefore  arc  not  representative  of  the  currents  descrip¬ 
tive  of  the  ODMOS.  For  example,  Station  C-33  shows  the  strong  ebb- flow  jet 
directed  at  an  angle  of  approximately  110  to  130  deg.  The  flood  flow  is  also 
shown  in  the  plot;  however,  the  direction  varies  somewhat  around  the  entrance 
of  the  Jetty.  Data  from  both  Stations  C-24  and  C-34  are  free  of  the  localized 
effects  of  the  jetty.  Both  vector  plots  indicate  a  tidal  ellipse  with  the 
major  axis  oriented  to  the  northeast- southwest  and  the  minor  axis  approxi- 
Buitely  90  deg  offset.  Since  these  two  gage  locations  can  be  considered  to  be 
representative  of  the  ODMDS,  the  uniformity  of  the  currents  with  respect  to 
depth  is  investigated  to  show  that  the  use  of  a  depth -averaged  velocity  in  the 
numerical  modeling  effort  is  justified.  Vector  plots  of  each  of  the  time 
series  for  each  of  the  depths  indicated  in  Table  2  for  Station  C-24  are  given 
in  Figure  12.  Results  show  similar  trends  at  each  depth  with  an  increase  in 
velocity  from  the  bottom  to  the  top. 

17.  Results  of  the  harmonic  analyses  as  well  as  the  continuous  vector 
diagrams  indicate  that  a  velocity  field  with  the  characteristics  of  those 
shown  for  gages  at  Stations  C-24  or  C-34,  i.e.  maximum  velocity  magnitudes  on 
the  order  of  40  to  50  cm/sec,  is  representative  of  current  patterns  immedi¬ 
ately  north  of  the  ODMDS.  In  both  cases,  residual  currents  are  directed  to 
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th*  northMUit.  Analysis  of  ths  onsite  neer*botton  nounted  EPA  gage  datemined 
whether  currents  within  the  OOKDS  are  significantly  different  frosi  those  at 
NQAA  sites  C-2A  and  C-34. 

18.  The  EPA  current  gage  is  located  in  the  northernmost  comer  of  the 
ODMDS,  as  shown  in  Figure  6.  The  instrument,  an  ENDECO  Type  174SSM  solid 
state  current  meter,  is  located  approxisMtely  6  ft  from  the  bottom.*  Data 
were  collected  during  the  periods  of  21  June -17  July,  24  July- 20  Aug,  and 

15  Nov- 11  Dec  1990.  Current  values  with  the  superimposed  low-frequency  com¬ 
ponent  and  the  directional  histograms  are  shown  in  Figures  13,  14,  and  15. 
Continuous  vector  plots  are  sho%m  in  Figure  16.  Results  of  the  harmonic 
analyses  of  Che  three  data  series  are  presented  in  Table  4. 

19.  Analysis  of  Che  data  presented  in  Table  4  and  Figures  13  throu^  16 
indicates  that  currents  at  the  EPA  gage  site  are  similar  in  both  magnitude  and 
orientation  to  those  of  NOAA  Stations  C-24  and  C-34.  For  example,  the 
coBd>ined  series  vector  plot  shown  in  Figure  16  is  similar  in  qualitative 
details  to  the  near -bottom  vector  plots  for  NOAA  Station  C-24,  shown  in  Fig¬ 
ures  12a  and  b.  It  is  therefore  reasonable  to  assume  that  the  depth-averaged 
harmonic  constituents  coiaputed  for  Station  C-24  can  be  used  to  characterize 
Che  tidal  forcing  at  the  Charleston  ODMDS. 

20.  The  harmonic  analyses  indicate  that  approximately  50  to  60  percent 
of  the  currents  can  be  attributed  to  tides.  To  account  for  the  source  of  the 
variance,  a  low-frequency  component,  similar  in  magnitude  to  those  shown  in 
Figures  8,  10,  and  12  through  14,  is  introduced  to  simulate  the  low-frequency 
component  of  the  observed  time  series. 

21.  The  long-term  modeling  goal  is  to  generate  a  data  base  of  simulated 
current  data  that  is  realistically  representative  of  currents  at  the  disposal 
site.  In  the  same  manner  that  the  wave  fields  were  simulated  to  reflect  the 
observed  statistical  distribution  of  the  WIS  hindcast  data,  a  time  series  of 
ctirrenCs  was  prepared  for  the  ODMDS,  based  on  the  harmonic  constituents  for 
NQAA  Station  C-24,  shown  in  Table  3.  Although  the  data  are  not  of  sufficient 
length  for  a  reliable  harmonic  analysis,  the  procedure  provides  an  approxiswte 
estimate  of  tidal  influence.  Inspection  of  the  low-  and  high-frequency 
portions  of  the  velocity  magnitude  as  well  as  the  actual  U  and  V  components  of 
the  data  shown  in  Figures  8  and  10  suggest  that  the  addition  of  a  long-period. 


*  Personal  Communication,  Dec  1990,  Philip  Murphey,  Physical  Scientist, 
Environmental  Protection  Agency. 
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c.  Series  3  d.  Series  1,2,  and  3 

Figure  16.  Current  vector  dlagrans 
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Table  4 


Harnontc  ATHIIYIU  nuta 


Series  #1  Series  #2 _  _ SgriffS  fl 


21  June 

-  17  July 

24  July  - 

20  Aug 

15  Nov  - 

11  Dec 

Const.* 

V 

Y 

u 

Y 

V 

Y 

Ma 

10.40 

6.25 

7.88 

5.44 

7.41 

5.21 

Sa 

1.19 

0.90 

1.43 

0.67 

0.77 

1.26 

Na 

2.70 

1.94 

1.44 

1.15 

1.69 

2.04 

Ki 

5.22 

1.87 

1.44 

3.24 

0.38 

0.44 

0i 

0.63 

0.61 

0.61 

0.83 

0.07 

0.56 

Pi 

4.17 

3.26 

1.31 

2.23 

0.17 

0.50 

RV 

0.55 

0.30 

0.50 

0.26 

0.68 

0.48 

Ave. 

3.21 

7.75 

-2.10 

0.33 

-1.11 

1.29 

*  Const.  -  constituent 
Ave.  -  average 

large>ainplitude  component  to  the  tidal  signal  would  produce  fluctuations  in 
the  simulated  current  time  series  that  would  be  representative  of  prototype 
conditions.  Therefore,  a  synthetic  U  and  V  tidal  component  with  an  amplitude 
of  20  cm/sec  and  a  period  of  8  days  was  added  to  the  constituent  list  shown  in 
Table  3.  The  resulting  tidal  signal  and  vector  is  shown  in  Figure  17.  Note 
that  the  magnitude  envelopes  are  similar  in  structure  to  the  observed  proto¬ 
type  data.  Therefore,  the  tidal  constituents  listed  in  Table  3  and  the  8-day 
low-frequency  component  are  used  to  simulate  tidal  hei^t  and  current  fluctua¬ 
tion  in  the  long-term  modeling  effort.  A  residual  current  of  U  -  3.9  and 
V  «  2.6  cm/sec  was  included  in  the  synthesized  tidal  signal  (values  calculated 
for  C-24). 

22.  A  single-value  velocity  is  specified  for  the  short-term  modeling 
effort  since  the  model  simulations  correspond  only  to  a  single  disposal  To 
represent  "worst -case  scenarios,"  maximum  envelope  current  values  were 
selected  to  simulate  the  maximum  spatial  excursion  of  suspended  sediment  fol¬ 
lowing  release  from  the  barge.  In  view  of  the  magnitudes  shown  in  Figure  11, 
two  depth -averaged  current  values  were  specified  for  the  short-term  simula¬ 
tions.  A  value  of  45.7  cm/sec  was  specified  for  the  northeast- southwest 
oriented  major  axis,  and  a  value  of  30.5  cm/sec  was  specified  for  the 
northwest-southeast  oriented  minor  axis  of  the  current  ellipse. 
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PART  III:  SHORT-TERM  MODELING 


General 

23.  The  short- terai  nodellng  component  of  this  investigation  examines 
the  immediate  effect  of  the  actual  disposal  operation  on  the  surrounding  area. 
Numerical  simulations  of  the  discharge  are  used  to  determine  whether  the 
combined  effects  of  the  local  topography  at  the  site  and  the  depth -averaged 
velocity  field  pose  a  threat  to  the  effectiveness  of  the  dredged  material 
disposal  operation.  Can  the  material  be  physically  placed  within  the  limits 
of  the  designated  site  as  the  material  descends  throu^  the  water  column  to 
the  ocean  floor,  or  are  the  local  currents  of  sufficient  magnitude  to  trans¬ 
port  material  out  of  the  site  before  deposition? 

24.  The  short-term  site  evaluation  phase  is  made  by  numerically 
modeling  the  disposal  operation  using  the  DIFID  numerical  model.  Tlieory  and 
background  of  the  model  are  reported  by  Johnson  and  Holiday  (1978) ,  Johnson 
(1990),  and  Johnson,  Trawle,  and  Adamec  (1988).  Applications  of  the  model  are 
reported  by  Trawle  and  Johnson  (1986) ,  Scheffner  (1992  and  in  preparation) , 
and  Schexfuer  and  Swain  (In  preparation) .  The  model  computes  the  time  history 
of  a  single  disposal  operation  from  the  time  the  dredged  material  is  released 
from  the  barge  until  it  reaches  equilibrium  on  the  ocean  floor.  The  DIFID 
model  separates  the  dumping  operation  into  three  distinct  phases.  In  the 
first  phase,  material  released  from  the  bln  is  assumed  to  form  a  hemlspheri- 
cally  shaped  cloud,  which  descends  throu^  the  water  column  under  the  influ¬ 
ence  of  gravity.  This  phase  is  called  the  convective  descent  phase. 

25.  The  convective  descent  phase  continues  until  the  cloud  of  material 
impacts  the  bottom  or  reaches  a  stable  point  of  neutral  buoyancy.  In  either 
case,  horizontal  spreading  of  material  marks  the  beginning  of  the  dynamic 
collapse  phase  in  which  the  material  spreads  horizontally.  When  the  rate  of 
spreading  becomes  less  than  spreading  due  to  turbulent  diffusion,  the  final 
transport -diffusion  phase  of  transport  begins.  An  idealization  of  all  three 
phases  of  the  short-term  disposal  are  shown  in  Figure  18. 
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26.  The  DIFID  model  requires  site-specific  input  data  to  quantitatively 
predict  the  short-term  sediment  fate  of  a  disposal  operation.  These  data 
include  the  physical  dimensions  of  the  dredge,  a  description  of  the  local 
environment  (depth  and  velocity  field) ,  a  knowledge  of  the  composition  and 
characteristics  of  the  material  in  the  dredge,  and  specification  of  the 
numerous  modeling  parameters  and  coefficients. 

27.  In  this  particular  case,  inner  harbor  dredged  material  was  sepa¬ 
rated  into  two  groups  based  on  sediment  size.  The  coarse  and  fine  material 
had  a  constitutent  breakdown  of  15/85  and  85/15  percent  silt- clay/  sand, 
respectively.  Additionally,  the  two  main  components  of  the  mean  depth- 
averaged  current  system  (northeast- southwest  and  northwest-southeast)  were 


applied  to  the  model.  Based  on  the  combination  of  various  model  input 

parameters,  four  case  studies  are  performed.  They  include: 

Case  1:  Worst-case  depth -averaged  northeast -southwest  current 
speed  of  1.5  ft/sec  and  fine-grained  sediment. 

b.  Case  2:  Worst-case  depth- averaged  northwest -southeast  current 
speed  of  1.0  ft/sec  and  fine-grained  sediment. 

£.  Case  3:  Worst-case  depth- averaged  northeast -southwest  current 
speed  of  1.5  ft/sec  and  coarse-grained  sediment. 

Case  4:  Worst-case  depth- averaged  northwest- southeast  current 
speed  of  1.0  ft/sec  and  coarse-grained  sediment. 

28.  Model  input  requires  specification  of  the  size  and  capacity  of 
the  dredge.  It  is  anticipated  that  70  percent  of  dredging  operations  will 
involve  a  "Manhattan"  class  dredge,  or  one  of  similar  dimensions.*  This 
dredge  will  deposit  material  in  approximately  45  to  55  ft  of  water.  Because 
the  DIFID  model  assumes  a  single  load  disposal,  it  is  assumed  that  the  entire 
volume  of  a  Manhattan  class  hopper  dredge  is  released  at  once.  Althou^  this 
assumption  may  not  be  entirely  accurate,  results  give  a  worst  case  scenario 
from  which  disposal  decisions  can  be  based.  Capacities  and  dimensions  of  the 
Manhattan  class  dredge  are  given  in  Table  5. 


Table  5 

CaPflCitlgs  and  Dimenslaps  of  the  Manhattan  Class  Dredge 


Overall  length  280  ft 

Width  52  ft 

Depth  21.5  ft 

Unloaded  draft  8  ft 

Loaded  draft  of  vessel  19.4  ft 

Vol\ime  3,000  cu  yd 


29.  Additional  site-specific  parameters  include  the  specification 
of  grid  resolution  Ax  ,  total  simulation  duration  T^  ,  and  the  long-term 
time -step  At  .  Parameter  estimation  is  presently  based  on  a  review  of  pre¬ 
vious  case  studies  and  the  ongoing  DIFID  model  sensitivity  investigation. 
Values  for  the  internal  model  coefficients  were  based  on  recommendations  and 


*  Personal  Communication,  Dec  1990,  B.  Kyzer,  USAE  District,  Charleston, 
Charleston,  SC. 
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Table  6 


Model  Input  Parameters  and  Coefficients 


Variables _  _ 


Grid  size,  Ak  (ft) 

100 

Number  of  cells: 

Cross -shore  direction 

105 

Alongshore  direction 

28 

Time-step,  At  (sec) 

80 

Duration  of  slsiulatlon,  T^  (sec) 

8,000  (fine-grained  site) 

Duration  of  slsnilatlon,  T^  (sec) 

8,000  (coarse-grained  site) 

Ambient  velocity  (cm/sec) 

Cases  1  and  3 

45.7 

Cases  2  and  4 

30.5 

Local  depth  (ft) 

65.0  (fine-grained  site) 

X-Dlrection  (on-offshore) 
bottom  slope  (deg) 

65.0  (coarse-grained  site) 

0.00 

Y- Direct ion  (alongshore) 
bottom  Slope  (deg) 

0.00 

Aa^lent  density  (gm/cc) 

1.018 

DINCRl 

1.0 

DINCR2 

1.0 

Entrainment  coefficient  AIAPHO 

0.235 

BETA 

0.0 

CM 

1.0 

Drag  coefficient  for  sphere,  CD 

0.5 

GAMA 

0.25 

Drag  coefficient  for  elliptic 
cylinder,  CDRAG 

1.0 

CERIC 

0.01 

CD3 

0.10 

CD4 

1.00 

ALPHAC 

0.0010 

Bottom  friction,  FRICTN 

0.0100 

FI 

0.10 

ALAMDA 

0.005 

AKYO 

0.05 

applications  reported  by  Johnson  (1990)  and  Johnson  and  Holiday  (1978).  The 
parameters  and  coefficients  tised  in  both  simulations  are  shown  in  Table  6. 

30.  Final  input  to  the  DIFID  model  is  the  specification  of  the  composl 
tion  of  the  solid  material  in  the  dredge  according  to  percent  volume  of  sand, 
clay  and  silt,  clumps,  rocks,  etc.  Each  cosiponent  must  be  defined  according 
to  its  respective  density,  concentration  by  volume,  fall  velocity,  and  voids 
ratio.  The  percent  distribution  of  the  coarse  sediments  was  an  85/15 -percent 
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distribution  of  ssnd/silt-clsy,  whereas  the  fine  sediaents  contained  a 
15/64/21 -percent  distribution  of  sand/silt-clay/clay-cluBps.*  These  percent¬ 
ages  represent  only  the  solids  portion  of  the  material.  The  total  fluid 
composition  of  each  sample  was  based  on  a  separate  percent  distribution 
computation  for  the  water  content  of  the  sand  portion  and  the  silt-clay 
portion.  The  coarse  materials  were  defined  to  contain  68-percent  solids,  of 
idilch  85  percent  is  sand  and  15  percent  is  silt-clay.  The  fine-grained  sam¬ 
ples  were  defined  as  42 -percent  solid,  with  15 -percent  sand  and  85 -percent 
silt-clay.  The  silt-clay  component  is  fvtrther  separated  into  a  mixture  of  75- 
percent  suspendable  material  and  25 -percent  clumps.  Final  results  of  the 
coaq>utatlons  are  sho%m  in  Table  7  for  the  fine-grained  material  and  in  Table  8 
for  the  coarse-grained  material. 


Table  7 

Fine-Grained  Sediment  Composition  and  Characteristics 


Density 

Volume 

Descrintion 

g/cc 

percent 

Sand 

2.600 

6.3 

Silt-Clay 

2.600 

26.8 

Clumps 

2.600 

8.9 

Water 

1.018 

58.0 

Table  8 

Coarse-Grained  Sediment  Composition  and  Characteristics 


Density 

Volume 

g/cc 

percent 

Sand 

2.600 

57.8 

Silt -Clay 

2.600 

10.2 

Water 

1.018 

32.0 

*  Personal  Communication,  Dec  1990,  B.  Kyzer,  USAE  District,  Charleston, 
Charleston,  SC. 
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praaantad  In  tha  following  sactlons. 
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Sbartrltni  Modftl  Siauldtlona 


31.  The  short-term  model  Is  esiployed  in  the  selection  of  a  disposal 
point  within  the  boundaries  of  the  designated  disposal  site  that  would 
minimize  any  disposal -related  effect  to  the  hard  bottom  areas.  For  this  pur¬ 
pose,  two  properties  of  the  disposal  operation  are  investigated:  (a)  the  con¬ 
centration  and  transport  characteristics  of  the  sediment  plume,  and  (b)  the 
maximum  horizontal  distance  from  the  disposal  point  at  which  an  appreciable 
aaraunt  of  material  settles,  i.e.  the  measureable  limit  of  deposition.  These 
two  properties  of  the  impact  are  addressed  by  the  model  during  both  the 
dredged  material  descent  and  collapse  phases.  During  the  descent  phase,  the 
size  of  the  sediment  cloud  increases  (diffuses)  and  consequently  becomes  less 
concentrated.  Calculations  during  this  phase  can  be  used  to  estimate  the  time 
change  in  sediment  concentration  with  depth  and  distance  from  the  barge.  Both 
concentration  distribution  and  total  deposition  results  are  presented 
separately  for  the  fine-  and  coarse-grained  disposal  operations.  Lastly,  due 
to  DIFID  model  constraints,  a  disposal  depth  of  65  ft  is  employed  in  the  simu¬ 
lation.  By  modeling  the  disposal  operation  in  a  65-ft  depth  rather  than  the 
measured  45-  to  55- ft  depth,  a  conservative  suspended  sediment  concentration 
estisuite  is  effectively  ensured. 


Fine-Grained  Disposal  Site  Analysis 

32.  To  simulate  the  fine-grained  disposal  operation  in  the  two  worst- 
case  current  conditions,  the  coefficients  presented  previously  were  input  to 
the  numerical  model  for  each  of  the  major  and  minor  current  axes.  Model 
results  include  the  spatial  distribution  of  the  sand  and  silt -clay  components 
of  the  sediment  load  in  the  form  of  a  nondimens ional  ratio  of  volume  of  solids 
to  volume  of  solution.  To  convert  these  values  to  mg/i,  multiply  by  sediment 
density  (2.60)  x  10*.  The  clay-clusq)  component  of  the  fine  material  settled 
to  the  sea  bottom  relatively  quickly;  consequently,  a  concentration  distribu¬ 
tion  is  not  provided.  An  example  of  transport  and  diffusion  of  the  sediment 
cloud  is  shown  in  Figures  19  through  21,  in  which  the  horizontal  distribution 
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of  the  svMpended  sedinent  concentration  of  the  silt-clay  cloud  is  shown  at  the 
30- ft  depth  (below  the  surface)  for  the  postduap  times  of  2,000,  4,000,  and 
6,000  sec.  These  concentration  snapshots  show  the  increase  in  size  and 
corresponding  decrease  in  concentration  of  the  settling  cloud  as  it  is  dis¬ 
persed  and  diffused  from  the  point  of  release  from  the  barge.  The  point  of 
disposal  in  Figures  19-21  is  at  grid  cell  10  (0.19-mile  point). 

33.  Results  of  the  concentration  computation  are  used  to  produce  a 
volume  rat:o  concentration  above  ambient  conditions  versus  distance  relation¬ 
ship  along  the  central  axis  of  the  grid  at  three  discrete  depths  for  three 
specified  time  periods  (i.e.,  along  the  axis  of  symmetry  at  grid  14  of  Fig¬ 
ures  19-21).  Quarter-points  were  selected  to  show  results  at  one-fourth,  one- 
half,  and  three-fourths  of  the  termination  time  following  the  initial  release 
of  material  from  the  barge.  These  plots  are  presented  only  for  the  silt-clay 
component  of  the  disposed  material  since  the  concentration  of  suspended  sand 
was  comparatively  insignificant.  In  Figures  22  and  23,  the  concentration 
history  plots  for  the  silt-clay  component  under  the  influence  of  the  worst- 
case  Man  current  speeds  of  1.5  (Case  1)  and  1.0  (Case  2)  ft/sec  are  pre¬ 
sented  ,  respectively . 

34.  Results  shown  in  Figures  22  and  23  represent  time-concentration 
histories  along  the  suspended  sediment  cloud  axis.  The  three  concentration 
profiles  shown  at  the  30- ft  level  of  Figure  22  correspond  to  the  central  axis 
of  Figures  19  through  21.  The  three  depths  of  30,  40,  and  50  ft  were  used  to 
demonstrate  the  sediment  distribution  through  the  water  column.  For  example, 
simulations  of  the  disposal  operation  indicate  that  very  little  suspended 
sediment  exists  in  the  upper  40  ft  of  the  water  column  6,000  sec  after  the 
initial  dump;  i.e.,  the  material  has  passed  through  that  depth.  The  examples 
presented  in  Figures  22  and  23  indicate  that  the  maximum  silt-clay  concentra¬ 
tion  is  located  in  the  middepth  region  and  that  the  silt-clay  concentration 
decreases  with  both  time  after  disposal  and  distance  from  the  release  point. 

A  summary  of  the  Cases  1  and  Case  2  silt-clay  concentration  simulations  are 
shown  in  Tables  9  and  10,  respectively. 

35.  In  Figures  24  and  25,  the  total  sediment  deposition  versus  distance 
along  the  axis  of  the  disposal  grid  is  shown  for  Cases  1  and  2,  respectively. 
Additionally,  a  three-dimensional  view  of  the  Cases  1  and  2  deposition  pattern 
is  shown  in  Figures  26  and  27,  respectively,  with  the  corresponding  contour 
plots  shown  in  Figures  28  and  29. 
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Plgura  19.  Suapandad  aadlaant  cloud  at  30-ft  depth, 
2,000  aec  after  diaposal 
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Figvre  20.  Suspended  sediaent  cloud  at  30* ft  depth, 
4,000  sec  after  disposal 
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Figure  22.  Silt-clay  concentration 
evolution  relationships  for  the 
fine-grained  site  (Case  1) 
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Figure  23.  Silt-clay  concentration 
evolution  relationships  for  the 
fine-grained  site  (Case  2) 


Table  9 

Simwjirv  of  Computed  MaxlBum  Suspended  Silt-Clav  Concentration.*  Case  I 


Time .  sec/Approximate  Distance 

from  Disposal,  miles 

Depth,  ft 

^99Q/Q,?? 

4.000/1.13 

6.000/1.75 

30 

5.8  X  10*5 

1.2  X  10'5 

4.0  X  10-® 

40 

5.3  X  10-5 

1.1  X  10-5 

3.6  X  10'® 

50 

3.3  X  10‘* 

7.1  X  10-® 

2.3  X  10'® 

*  Concentration 

in  volume  ratio  of  solids 

to  solution  i 

above  ambient. 

Table  10 

SiiminuT-Y  9f  Computed  Maximum  Suspended  Silt-Clav  Concentration.*  Case  2 


Time,  sec /Approximate  Distance  from  Disposal,  miles 


2. 000 /p, 37 

4.000/0.75 

6.000/1.12 

30 

4.7  X  10-5 

1.0  X  10-5 

3.5  X  10'® 

40 

4.2  X  10-5 

1.0  X  10-5 

3.2  X  10-* 

50 

2.7  X  10-5 

6.0  X  10'® 

2.0  X  10’® 

*  Concentration  in  volume  ratio  of  solids  to  solution  above  ambient. 


DISTANCE  IN  MILES 

Figure  24.  Total  deposition  pattern 
for  the  fine-grained  site  (Case  1) 
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Flgurs  25.  Total  dapoaition  pattern 
for  the  fine-grained  site  (Case  2) 


Figure  26.  Three -dioens tonal  vlev 
of  the  fine-grained  site  deposition 
pattern  (Case  1) 
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Figure  27.  Three -diaensional  view 
of  the  fine-grained  site  deposition 
pattern  (Case  2) 
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Figure  28.  Contour  plot  of  the  fine-grained 
site  deposition  pattern  (Case  1) 
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Figure  29.  Contour  plot  of  the  fine 'grained 
site  deposition  pattern  (Case  2) 

36.  To  aid  the  engineer  or  scientist  in  selecting  a  disposal  location 
that  would  minimize  the  effect  on  hard  bottom  areas,  pliime  concentration 
information  has  been  compiled  and  is  presented  in  Figure  30.  This  figure 
describes  the  evolution  of  depth 'averaged  peak  plume  concentration  as  it  is 
advected  from  the  disposal  point.  In  the  figure,  exponential -type  regression 
curves  are  given  for  Cases  1  and  2;  the  corresponding  equation  follows. 

•*<1  or  2) 


Where  Q  is  the  allowable  silt-clay  s\ispended  sediment  volume  ratio  concen¬ 
tration  X  10^,  and  a  and  ^  are  equal  to  -2.2  and  116  for  Case  1,  and  -2.8 
and  71  for  Case  2. 

37.  A  simple  method  that  can  be  used  to  determine  whether  or  not  a 
particular  disposal  location  meets  the  required  minimal  hard  bottom  effect 
requirements  is  given: 

£.  Determine  an  acceptable  suspended  sediment  concentration  level 
at  the  location  of  the  hard  bottom. 
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40 


Figure  30.  Evolution  of  depth -averaged  peak  plume  concen¬ 
tration  for  Cases  1  and  2  of  fine-grained  material 


Determine  the  value  of  Xx  and  X2  from  Figure  31  or  from 
Equation  3. 

£.  Construct  an  ellipse  about  the  selected  disposal  point  (center) 
with  a  major  axis  of  length  2Xx  running  in  the  northeast - 
southwest  direction  and  a  minor  axis  of  length  2X2  running  in 
the  northwest -southeast  direction. 

j|.  If  the  hard  bottom  areas  fall  outside  the  boundary  of  the 

ellipse,  the  selected  disposal  location  is  acceptable;  if  not, 
select  another  disposal  point  and  reapply  this  test. 


Coarse-Grained  Disposal  Site  Analysis 


38.  The  single-load  deposition  sinnilation  for  the  coarse-grained  mate¬ 
rial  was  performed  using  the  coefficients  shown  in  Tables  6  and  8.  Results  of 
the  sinnilations  showed  that  the  coarsest  fraction  of  this  material  descended 
rapidly  to  the  ocean  floor,  leaving  no  sand  in  sxispension  within  the  water 
coliimn.  The  silt-clay  fraction  of  this  material,  however,  rapidly  moved  into 
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•usp«nalon  as  can  be  seen  in  Figures  31  and  32.  As  in  the  case  of  the  fine* 
grained  uterial,  the  concentration  of  silt-clay  slightly  increases  with 
distance  above  the  bottom  and  decreases  with  time.  However,  the  siaxlmum  value 
of  coarse-grained  suspended  sediment  concentration  is  an  order  of  magnitiide 
less  than  that  of  the  fine-grained  material.  A  summary  of  the  Case  3  and 
Case  4  silt-clay  concentration  simulations  are  shotm  in  Tables  11  and  12, 
respectively.  In  Figures  33  and  34,  the  total  sediment  deposited  versus  dis¬ 
tance  along  the  axis  of  the  disposal  grid  is  shown  for  Cases  3  and  4, 
respectively. 

39.  Additionally,  a  three-dimensional  view  of  the  Cases  3  and  4  deposi¬ 
tion  pattern  is  shown  in  Figures  35  and  36,  respectively,  with  the  correspond¬ 
ing  contour  plot  shown  in  Figures  37  and  38.  Based  on  these  figures,  the 
maximum  thickness  of  deposition  for  Cases  3  and  4  is  approximately  0.45  and 
0.5  ft,  respectively.  The  deposition,  for  both  cases,  covers  an  approximate 
700-ft-diameter  area.  Deposition  is  shown  to  be  confined  to  this  immediate 
area. 

40.  As  in  the  case  of  the  fine-grained  material,  plume  concentration 
Information  for  the  coarse-grained  site  has  been  compiled  and  is  presented  in 
Figure  39.  This  figure  describes  the  evolution  of  depth -averaged  peak  plume 
concentration  as  it  is  advected  from  the  disposal  point.  In  the  figure 
exponential -type  regression  curves  are  given  for  Cases  3  and  4;  the  corre¬ 
sponding  equation  follows. 


X 


3  or  4 


(4) 


where  C  is  the  allowable  silt-clay  suspended  sediment  volume  ratio  concen¬ 
tration  X  10^,  and  a  and  ^  are  equal  to  -2.5  and  35  for  Case  3,  and  -4.8 
and  52  for  Case  4. 

41.  In  both  the  fine-  and  coarse-grained  DIFID  analyses,  depth -averaged 
velocities  of  1.5  and  1.0  ft/sec  were  assumed  as  a  conservative  estimate  of 
the  northeast-southwest  and  northwest- southeast  velocity  component,  respec¬ 
tively.  As  shown  in  the  prototype  data  analysis,  these  velocity  magnitudes 
represents  a  much  higher- than-average  condition,  and,  as  such,  the  results 
presented  for  the  short-term  simulation  may  be  considered  conservative  with 
respect  to  the  dispersion  of  the  suspended  sediments.  An  analysis  of  the 


53 


DISTANCE  IN  MLES 
c.  50- ft  depth 

Figure  32.  Silt-clay  concentration  evolution 
relationships  for  the  coarse-grained  site 

(Case  4) 


Table  11 


g.«.^rv  ftf  the  CoBDUtad  Maxlaua  Suananded  Sllt-Clav  Concentration.*  C4M  3 


Tine.  Dlatance  fron  Disposal,  nlles _ 

Depth,  ft 

2>flflQ/g.50 

4.MQ/1.12 

6.000/1.70 

30 

1.3  X  10-3 

2.8  X  10-* 

9.2  X  10-^ 

40 

1.2  X  10-3 

1.5  X  10-* 

8.4  X  10-’ 

50 

7.7  X  10-® 

1.6  X  10-* 

5.1  X  10-^ 

*  Concentration  in  voluae  ratio  of  solids  to  solution  above  a8A>lent. 


Table  12 

Sxnimarv  of  Computed  Maxlaua  Suspended  Sllt-Clav  Concentration.*  Case  4 


_ Tine.  sec/Approxlnate  Distance  fron  Disposal,  nlles _ 

Pvpth..  ft 

2.000/0.54 

4.  PM/0 .7^ 

^.PQQ/i.l4 

30 

1.0  X  10-3 

2.3  X  10-* 

7.8  X  10-’ 

40 

1.0  X  10-3 

2.1  X  10-* 

7.1  X  10-^ 

50 

6.3  X  10-® 

1.3  X  10-® 

4.5  X  10-^ 

*  Concentration  In  volune  ratio  of  solids  to  solution  above  aid>lent. 


short •tern  analysis  -results  will  be  presented  following  the  long- tern  slsnila- 
tlons  described  In  Part  IV. 
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l.SO  1.75 


Flgur*  36.  Thr»»>diMi»ioiial  vlav  of  tho 
coarse -grained  site  deposition  pattern 
(Case  4) 


Figure  37.  Contour  plot  of  the  coarse-grained 
site  deposition  pattern  (Case  3) 
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Figure  39.  Evolution  of  depth-averaged  peak  plune  concentration 
for  Cases  3  and  4  of  fine-grained  material 

where  C  is  the  allowable  silt-clay  suspended  sediment  volume  ratio  concen¬ 
tration  X  10®,  and  a  and  /3  are  equal  to  -2.5  and  35  for  Case  3,  and  -4.8 
and  52  for  Case  4. 

In  both  the  fine-  and  coarse-grained  DIFID  analyses,  depth-averaged 
velocities  of  1.5  and  1.0  ft/sec  were  assumed  as  a  conservative  estimate  of 
the  northeast -southwest  and  northwest -southeast  velocity  component,  respec¬ 
tively.  As  shown  in  the  prototype  data  analysis,  these  velocity  magnitudes 
represents  a  much  higher -than- average  condition,  and,  as  such,  the  results 
P^®®®ntad  for  the  short-term  simulation  may  be  considered  conservative  with 
respect  to  the  dispersion  of  the  suspended  sediments.  An  analysis  of  the 
short-term  analysis  results  will  be  presented  following  the  long-term  simula¬ 
tions  described  in  Part  IV. 
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PART  IV:  LONG-TERM  MOOEUMG 

General 

42.  The  long- tern  slnulatlon  phase  of  the  site  characterization  study 
investigates  the  behavior  of  a  dredged  naterial  mound  over  time.  This  analy¬ 
sis  is  accoaplished  by  developing  a  means  of  classifying  disposal  sites  as 
either  dispersive  or  nondispersive  based  on  %diether  local  wave  and  velocity 
fields  are  adeqxiate  to  erode  and  transport  significant  amounts  of  material 
from  the  site.  The  local  currents  can  be  due  to  normal  tidal  action  and  mean 
flow  circulation  patterns  or  storm- related  activity.  Sediment  transport 
calculations  tise  these  waves  and  currents  to  estimate  mound  stability  as  a 
function  of  the  local  bathymetry  and  sediment  characteristics  at  both  the 
fine-grained  and  coarse-grained  sites.  The  transport  relationships  used  for 
this  study  were  limited  to  noncohesive  sediments.  Although  a  majority  of  the 
fine  grained  sediments  are  classified  as  silt/clay,  they  contain  approximately 
15%  sand.  The  assumption  was  made  that  the  overall  material  behaves  in  a 
noncohesive  manner  and  can  be  approximated  with  noncohesive  relationships. 

This  assvunption  is  consistent  with  results  presented  by  Kamphius  (1990)  stat¬ 
ing  that  "If  sand  or  gravel  is  presented  in  the  eroding  stream  or  overlying 
the  cohesive  formation  in  a  discontinuous  layer,  the  design  should  be  based  on 
the  sediment  transport  characteristics  of  the  granular  material."  Details  of 
the  specified  grain  sizes  for  the  coarse  and  fine  grained  sediments  are  pre¬ 
sented  below  under  Input  Data  Requirements. 

43.  This  final  phase  of  the  site  evaltiation  represents  an  extension  of 
the  short-term  fate  analysis  of  Part  III  in  which  site  dispersiveness  was 
based  on  the  ability  to  effectively  place  material  within  a  designated  site 
during  the  disposal  operation.  The  long-term  analysis  begins  with  the  assump¬ 
tion  that  the  short-term  disposal  operation  is  successful  in  creating  a  stable 
mound  configuration.  Whether  the  mound  is  dispersive  or  nondispersive  depends 
on  whether  the  local  wave  and  current  conditions  are  capable  of  resuspending 
and  transporting  significant  amounts  of  material  from  the  mound  such  that 
areas  adjacent  to  the  disposal  site  are  affected. 

44.  The  long-term  site  stability  analysis  approach  adopted  for  this 
study  uses  the  simulated  wave  and  current  time  series  described  in  Part  II  to 
provide  a  quantitative  estimate  of  the  stability  of  the  mound  as  a  function  of 
localized  environmental  conditions.  The  analysis  approach  is  based  on  coupled 
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hydrodynamic  and  sediment  transport  models  that  compute  the  transport  of 
noncoheslve  sediment  as  a  function  of  the  local  velocity  and  depth.  The 
resulting  distribution  of  transport  is  used  in  a  sediment  continuity  model  to 
compute  changes  in  the  bathymetry  of  the  sediment  mound.  Bath3nBetry  change 
computations  are  made  at  every  3-hr  time- step.  The  long-term  simulations  of 
mound  stability  indicate  whether  the  local  wave  and  current  regime  at  the 
disposal  site  are  of  sufficient  siagnitude  to  suspend  and  transport  bottom 
sediments . 


Input  Data  Requirement 

45.  The  site  stability  methodology  is  dependent  on  the  accurate  predic¬ 
tion  of  sediment  transport  at  the  site  under  investigation.  Empirical 
relationships  for  computing  sediment  transport  as  a  primary  function  of  depth- 
averaged  water  velocity,  local  depth,  and  sediment  grain  size  were  reported  by 
Ackers  and  White  (1973).  A  modification  of  the  transport  equations  was  intro¬ 
duced  to  reflect  an  increase  in  sediment  transport  rate  when  ambient  currents 
are  accompanied  by  surface  waves.  The  modification,  in  the  form  of  an 
effective  increase  in  the  depth- averaged  current,  is  based  on  concepts 
developed  by  Bljker  (1967)  and  implemented  by  Swart  (1976) .  This  additional 
transport  reflects  the  fact  that  wave- induced  orbital  velocities  are  capable 
of  suspending  bottom  sediments,  independent  of  the  sediment  put  in  suspension 
by  mean  currents.  The  total  amount  of  sediment  put  into  suspension  by  waves 
and  currents  is  then  transported  by  the  ambient  current  field. 

46.  The  modified  Ackers -White  relationships  are  used  to  compute  the 
transport  of  uniformly  graded  noncoheslve  sediment  in  the  grain  diameter  (Dso 
for  example)  range  of  0.04  to  4.00  mm  (White  1972).  Consequently,  Djo  values 
of  0.100  and  0.063  mm  were  selected  to  represent  the  fine  sand  and  silt-size 
sediments  to  be  deposited  at  the  ODMOS.  Computed  sediment  transport  versus 
depth- averaged  velocity  for  a  range  of  depths  corresponding  to  those  at  the 
coarse-grained  site  are  shown  in  Figures  40  and  41  for  the  Dso  values  of  0.100 
and  0.063  mm,  respectively.  The  WIS  summary  mean  wave  height  of  0.5  m  and 
wave  period  of  3.9  sec  were  specified  in  the  preparation  of  this  family  of 
curves.  From  these  two  figures,  the  threshold  velocities  necessary  for  the 
initiation  of  sediment  erosion  is  almost  identical  for  both  sediment  sizes, 
l.e.  1.0<U<1.5  ft/sec. 
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Figure  40.  Sediment  transport 'velocity  relationships 

for  Djo  <-  0.10  BB 
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Figure  41.  Sediment  transport -velocity  relationships 
for  Dso  -  0.0625  on 
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47.  The  simulated  wave  and  current  boundary  condition  time  series  at 
the  disposal  site  are  used  to  evaluate  the  diapers ive/nondispersive  character¬ 
istics  of  the  site.  This  computation  requires  the  specification  of  the 
geometric  features  of  a  typical  or  existing  disposal  mound  located  within  the 
limits  of  the  designated  site.  The  disposal  procedures  followed  at  the 
Charleston  ODMDS  are  to  dispose  material  in  parallel  mounds  at  a  spacing  of 
approximately  300  ft.  Typical  mounds  have  a  width  of  approximately  1,000  ft 
and  a  length  of  up  to  5,000  ft.  Heights  of  the  elongated  features  are  limited 
so  that  depths  do  not  become  less  than  approximately  30  ft.*  Average  depths 
within  the  site  are  on  the  order  of  35  to  40  ft;  therefore,  a  local  depth  of 
37  ft  was  specified  for  the  long-term  simulations.  A  symmetric  design  mound 
configuration  of  1,000  by  1,000  ft  with  a  height  of  7  ft  was  specified  as  a 
means  of  predicting  maximum  erosion  rates.  For  example,  computed  migration 
rates  along  the  long  axis  of  the  elongated  mound  would  be  minimal  if  the  cur¬ 
rent  was  aligned  with  the  long  axis.  By  specifying  a  symmetric  mound,  maximum 
migration  rates  are  computed  because  directional  effects  are  minimized.  A 
three-dimensional  perspective  view  and  contour  map  of  the  test  mound  are  shown 
in  Figures  42  and  43. 


Long-Term  Model  Simulations 

48.  The  long-term  analysis  described  in  the  following  section  uses  wave 
and  velocity  time  series  to  compute  the  time  evolution  of  a  dredged  material 
mound.  A  quantitative  assessment  of  mound  stability  is  made  by  computing  the 
location  of  the  centroid  of  the  mound  along  the  central  mound  axis  for  each 
computational  time  step  of  the  simulation.  These  computations  are  made  by 
balancing  the  summation  of  moments  at  each  computational  grid.  Simulation 
results  are  also  presented  in  the  form  of  postsimulation  perspective  and  con¬ 
tour  plots  as  well  as  time  evolution  plots  of  the  changing  cross-sectional 
profile  along  the  axis  of  the  mound. 

49.  The  stability  analysis  is  made  by  estimating  mound  response  to  long 
periods  of  exposure  to  the  wave  and  current  conditions  described  in  Part  II. 

In  addition  to  this  normal  condition  simulation,  a  storm  event  analysis  was 
performed  in  an  attempt  to  investigate  single -event -related  erosion  of  the 


*  Personal  Communication,  Nov  1991,  S.  Morrison,  USAE  District,  Charleston, 
Charleston,  SC. 
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OPPSHORC 

Figure  43.  Idealized  disposal  Bound  contour  map 
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test  mound.  A  design  storm  event  was  selected  with  a  half-sine  wave  shape, 
with  a  duration  of  24  hr  and  a  peak  depth -averaged  velocity  of  2.0  ft/sec.  A 
storm  of  this  sustained  magnitude  and  duration  could  be  descriptive  of  a  mod¬ 
erate  northeaster  or  hurricane  event. 

50.  Two  long-term  simulations  were  made  to  investigate  the  stability  of 
a  disposal  mound  composed  of  two  different  sediment  sizes.  The  first  simula¬ 
tion  specified  a  grain  size  of  0.100  mm,  corresponding  to  a  fine  sand;  the 
second  specified  a  diameter  of  0.063  mm,  corresponding  to  a  noncohesive 
sediment  approaching  the  silt/clay  size.  A  2-mo  simulation  of  each  sediment 
size  was  performed  to  test  the  design  mound  stability  to  the  wave  field  shown 
in  Figure  3  and  the  current  fields  described  by  the  constituents  listed  in 
Table  3.  Results  of  both  simulations  indicate  that  the  moxind  experiences 
erosion  and  migration  under  the  forces  associated  with  normal  tidal  and  circu¬ 
lation  action.  A  migration  rate  of  approximately  15  and  60  ft/mo  were 
computed  for  the  0.100-  and  0.063-mffl  sediments,  respectively.  These  results 
show  that  the  long-term  direction  of  migration  is  to  the  northeast,  along  the 
major  tidal  current  ellipse  and  in  the  direction  consistent  with  computed 
average  flows.  Although  the  current  records  used  in  the  reconstruction  of  the 
current  time  series  were  based  on  data  of  limited  duration,  the  computed  mag¬ 
nitudes  and  directions  are  realistically  descriptive  of  normal  conditions  at 
the  ODMDS.  It  is  unlikely  that  longer  duration  current  data  would  show 
results  that  are  substantially  different  from  those  shown  in  the  present  anal¬ 
ysis.  Also,  the  computed  magnitude  and  direction  patterns  for  the  NOAA  sites 
are  consistent  with  comparable  computations  for  the  five  samples  at  the  EPA 
site.  Therefore,  it  can  be  reasonably  concluded  that  the  currents  shown  in 
Figure  14  are  representative  of  the  site  and  that  these  forcings  are  adequate 
to  transport  material  from  the  disposal  mound.  In  view  of  the  repetitive 
nature  of  the  velocity  field  and  the  fact  that  the  imposed  wave  field  corre¬ 
sponds  to  the  high  energy  winter  period  beginning  1  January  of  the  simulated 
year,  simulations  of  longer  than  2  months  were  not  necessary. 

51.  The  long-teiHB  simulation  indicates  that  the  test  mound  is 
dispersive  whenever  the  current  magnitudes  exceed  approximately  1.5  ft/sec; 
however,  normal  conditions  are  not  sufficiently  severe  to  cause  massive 
erosion.  Therefore,  the  mound  was  stibjected  to  velocity  and  wave  field  condi¬ 
tions  to  determine  the  levels  of  severity  of  erosion  that  could  be  anticipated 
from  a  storm  event  of  moderate  intensity.  The  specified  storm  could  be  a 
northeaster  or  hurricane;  therefore,  the  current  hydrograph  could  be  from  any 
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direction.  Rsaults  of  tho  scon  surge  slaulstion  show  the  sound  Co  ■igrsce 
155  fC,  ssny  Cises  che  erosion  rsCe  sssociaCsd  wlch  nonsal  waves  and  currenCs. 
The  fine-grained  sound  (0.063  ns)  was  used  for  Che  scon  sinulacions  Co 
represenC  a  worse-case  erosion  scenario.  A  ploC  of  Che  poscscon  sisulaCion 
perspecCive  and  conCour  map  of  Che  mound  and  che  compuced  cross -secCional 
evoluCion  of  Che  mound  axis  are  shown  in  Figures  44  and  45. 
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Figure  44.  SCon-evenC  simulacion  perspecCive  and  conCour  map 
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Figure  45.  Storm-event  simulation  mound  axis  migration  history 


PART  V:  CONCLUSIOMS 


52.  A  dispersion  analysis  of  the  Charleston  ODMDS  was  conducted  by 
CERC,  WES,  to  quantify  the  dispersive  character  associated  with  disposal  of 
dredged  material  at  the  existing  site.  Of  specific  concern  was  whether  or  not 
the  disposal  operation  could  pose  a  threat  to  recently  discovered  hard  bottom 
reef  communities  located  within  the  boundaries  of  the  existing  site.  The 
dispersion  analysis  was  conducted  in  two  phases.  First,  a  short-term  analysis 
was  performed  in  tdiich  the  immediate  effects  of  the  disposal  operation  were 
investigated  by  computing  the  spatial  limits  of  sediment  deposition  on  the 
ocean  bottom  and  the  time  rate  of  change  of  suspended  sediment  concentrations 
within  the  water  column.  Second,  a  long-term  analysis  was  conducted  to  deter¬ 
mine  whether  local  waves  and  currents  were  sufficiently  severe  to  erode  and 
transport  significant  amounts  of  material  from  the  existing  disposal  mound 
onto  the  hard  bottom  reef  areas.  Both  short-  and  long-term  analyses  are  based 
on  numerical  model  predictions  xising  site-specific  wave  and  current  boundary 
conditions  developed  throu^  analysis  of  existing  data.  The  phases  of  data 
analysis  and  short-  and  long-term  modeling  are  summarized  in  the  following 
section. 


Prototype  Data  Analysis 

53.  Although  the  Charleston  ODMDS  is  reasonably  dispersive,  the  goal  of 
this  investigation  was  to  quantify  the  direction  and  magnitude  of  sediment 
transport  with  respect  to  the  point  of  disposal.  This  goal  requires  an  inves¬ 
tigation  and  analysis  of  the  local  waves  and  currents  that  are  responsible  for 
eroding  and  transporting  material  from  the  site.  Wave  field  data  for  this 
purpose  were  obtained  from  the  Phase  III  WIS  20-year  hindcast  data  base  for  a 
WIS  station  located  offshore  of  the  entrance  to  Charleston  Harbor.  Two 
sources  of  velocity  data  were  located  that  were  of  sufficient  quality  that 
reliable  estimates  of  current  magnitude  and  direction  at  the  site  could  be 
computed.  The  sources  of  data  were  the  NOAA  and  the  EPA.  Both  data  sets  were 
analyzed,  and  comparable  results  were  obtained  indicating  that  currents  at  the 
site  are  tidally  dominated  with  a  maximum  major  axis  current  amplitude  of  40 
to  50  cm/sec  to  the  northeast- southwest  and  a  minor  amplitude  axis  of  approxi¬ 
mately  25  cm/sec  oriented  in  a  northwest-southeast  direction.  Mean  values 
were  on  the  order  of  less  than  10  cm/sec  in  a  generally  northeasterly 
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direction.  Although  the  data  used  for  this  analysis  were  limited  in  duration, 
the  results  are  sufficiently  accurate  for  the  dispersion  computations.  Since 
the  disposal  site  is  located  in  relatively  shallow  water,  i.e.  less  than  45  to 
55  ft,  the  occurrence  of  even  moderate  storm  events  can  have  a  greater 
erosional  impact  on  the  BK>und  than  those  attributed  to  normal  conditions.  For 
this  reason,  a  storm  event  analysis  was  included  to  quantify  the  degree  of 
long-term  storm- related  erosion. 

Short-Term  Analysis 

54.  Two  sediment  compositions  were  specified  in  the  short-term  disper¬ 
sion  analysis,  a  disposal  load  composed  of  predominantly  sand  and  one  of  pre¬ 
dominately  silt  and  clay.  Results  of  the  short-term  analysis  indicate  that 
both  the  sand  and  sllt/clay  materials  fall  rapidly  to  the  ocean  floor  and  do 
not  affect  regions  beyond  1.5  miles  from  the  point  of  disposal.  Computations 
of  the  movement  and  dispersion  of  the  suspended  sediment  cloud  at  the  1.5  mile 
point  indicate  that  maximum  depth-averaged  sediment  concentrations  in  ratios 
of  volume  of  solids  to  voltime  of  solution  above  ambient  (background)  condi¬ 
tions  are  on  the  order  of  3.0  x  10'*.  Concentrations  at  this  level  are  equiv¬ 
alent  to  7.8  mg/i  and  correspond  to  the  predominantly  sllt/clay  sediment.  The 
corresponding  value  for  the  fine  sand  is  7.5  x  10'^  (2.0  mg/i),  approximately 
1.5  miles  from  the  point  of  release.  Figures  30  and  39  provide  guidance  for 
locating  disposal  positions  that  minimize  any  suspended  sediment  impact  to  the 
live  reef  area.  As  shown  in  both  figures,  disposal  locations  can  be  located 
that  do  not  significantly  affect  the  reef  area.  For  example,  near-reef  dis¬ 
posal  could  be  made  during  ebb  flow  when  currents  are  to  the  northeast.  When 
currents  are  reef -directed,  neither  the  sand  nor  the  silt/clay  materials  pose 
a  problem  if  the  disposal  site  is  greater  than  approximately  1.5  miles  from 
the  reefs. 


Long-Term  Analysis 

55.  I«ng-term  simulation  of  disposal  mound  stability  was  based  on 
numerical  mound  stability  computations  based  on  site -specific  wave  and  current 
botindary  condition  data.  Two  design  sediment  mounds  were  tested,  one  cor¬ 
responding  to  a  fine  sand  (0.100  mm)  and  one  to  a  finer-grained  sediment 
(0.063  mm).  The  calculations  are  based  on  the  modified  Ackers-White  eq\iations 
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of  traaaport  for  mmeohoai^  sodlaont.  Siaulaciono  ahowod  that  both  aadiaant 
■ouada  ifora  dlaparaivo  with  raapact  to  noroal  wava  and  tidal/circulation  cur- 
ranta.  with  nigration  ratas  on  tha  ordar  of  15  and  60  faat/nonth  for  tha  aand 
and  fina>grainad  sitas,  raspactivaly.  Both  rasults  indicata  movanant  to  tha 
northaast,  consistant  with  conputad  currant  avaragas. 

56.  Sinca  tha  ODMDS  is  locatad  in  shallow  watar,  effacts  of  a  stom- 
valocity  fiald  wara  sinulatad  to  daaonstrata  tha  af facts  of  storm* related 
erosion.  Tha  sinnilation  of  a  aodarata- intensity  event  (2  ft/sec)  with  a  24-hr 
storm  showed  the  migration  of  a  0.063 -mm  noncohasive  sediment  mound  to  be 
approximately  155  ft.  Since  storm  currents  can  be  directed  in  virtually  any 
direction,  it  can  be  seen  that  long-term  erosion  can  easily  be  storm  domi¬ 
nated.  Coaiputed  migration  rates  do  not,  however,  indicate  rapid  and  massive 
erosion  that  would  affect  areas  far  removed  from  the  mound.  Therefore,  as 
long  as  the  disposal  mounds  are  located  a  reasonable  distance  to  the  northeast 
of  the  reef  area,  no  significant  long-term  effect-  to  the  reef  area  should  be 
experienced. 


57.  In  conclusion,  both  short-  and  long-term  analyses  of  the  Charleston 
ODMDS  shoved  that  the  site  is  moderately  dispersive  and  that  the  primary 
direction  of  dispersion  and  erosion  is  to  the  northeast.  Computations 
presented  herein  provide  guidance  that  should  aid  in  locating  specific  points 
of  disposal  of  material  which  would  minimize  the  threat  to  the  live  bottom 
habitat  area. 
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